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Abstract 

New Raman and incoherent neutron scattering data at various temperatures 
and molecular dynamic simulations in amorphous silica, are compared to ob- 
tain the Raman coupling coefficient C{uj) and, in particular, its low frequency 
limit. This study indicates that in the uj — > limit C(w) extrapolates to 
a non vanishing value, giving important indications on the characteristics of 
the vibrational modes in disordered materials; in particular our results indi- 
cate that even in the limit of very long wavelength the local disorder implies 
non-regular local atomic displacements. 
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The vibrational dynamics in topologically disordered systems is one of the most intriguing 
problems of present-day condensed matter physics In this respect, an important issue is 
the excess of vibrational states, that has been shown to exist in glasses at low frequencies 
by: (i) the specific heat that, above a few degrees Kelvin, exhibits a bump in Cp/T 3 
and (ii) the total density of states g(uj) that, at corresponding frequencies, shows a 
broad excess band in the plot of g{uo)/uo 2 ||-f| often referred to as boson peak (BP). Both 
findings indicate a non-Debye behaviour of the vibrational properties at frequencies below 
10-^50 cm" 1 , depending on the specific glass. The nature of the excess modes is still object 
of speculations [|TJ. In particular, in the case of S1O2 there are two prevailing hypotheses. 
According to the first one, strong phonon scattering by the structural disorder induces 
a localization of the vibrational states, and the excess of states appears at the crossover 
between the low frequency propagating and the high frequency localized modes ||. In the 
other picture, collective propagating modes persist up to frequencies higher than the BP 
frequencies, and the BP itself reflects their density of states PH|-|I2|. 

Spectral information on the BP is usually obtained by Inelastic Neutron Scattering (INS) 
or by Raman scattering (RS) spectra. In both cases, the first order scattering intensity, 
I(r,n)(u,T), is connected to the vibrational density of states by |I3fl : 



I(r,n)(u,T) cx n + 1 g(u)C( R , N )(uj)] (1) 

UJ 

here n{u) is the Bose population factor and Cm^{Lo) is the probe-excitations coupling 
function. The incoherent neutron scattering depends only on the absolute motion of atoms in 
space, so Cn(u) = 1. On the contrary, RS detects the relative displacement of neighbouring 
atoms and Cr(lo) turns out to be a complicated function of uj [fL3| . Since Cr(uj) is usually 



unknown, it results from Eq. 1 that RS spectra are not sufficient to determine g(uj), though 
this experimental technique is excellent for luminosity and resolving power. 

On a theoretical ground, different models have been developed for the frequency depen- 
dence of C(u). In particular, it was shown that i) for slightly distorted plane wave vibrations, 
C(u) oc uj 2 jnj; ii) within the framework of the soft potential model ||15|| , C{u) = const; and 



Ill) 



) the fracton-like model implies C(u) oc u a [ I6|,|17| . It is worth to note that all the models 



proposed so far, predict a power-law behaviour for C(u). On the contrary, the results of 
numerical simulations are much more complex and rarely produce scaling behaviour, even 



on simple model systems ||T8|-21 



On the experimental side, the direct comparison between RS and INS spectra appears 
to be a reliable procedure to determine C (u) , and then to check the validity of the proposed 
theoretical models. Though some attempts in this direction were made in the past a 



definite conclusion as to the shape of C(u) has not been reached |8 , 17| , 23|] , mainly due to the 



presence of "spurious" effects in RS (namely, the presence of quasi-elastic scattering (QES) 
or of luminescence). 

The aim of the present work is to determine the spectral shape of C (u) in vitreous silica 
in the acoustic frequency range (uo < 150 cm -1 ), at different temperatures from new RS and 
INS measurements. The luminescence is carefully analysed and eliminated from RS spectra. 
The Ir/In ratio results to be T-independent, except in the low frequency region where the 
effects of QES in RS are also evident. The coupling function, as defined in Eq.[TJfor harmonic 
excitations, is then extrapolated as the T — > limit of Ir/In- The obtained function does 
not show the power-law behaviour previously proposed in the theoretical models. Further, it 
appears not to vanish in the u — > limit, giving important information on the characteristics 
of the normal modes of vibration in disordered systems. These findings are also compared 
with the outcoming of a molecular dynamics study of vitreous silica, which confirms the 
experimental result and allows a direct inspection of the vibrational modes. 

The S1O2 suprasil sample purchased from Goodfellow was used for light scattering, while 
Heralux quartz from Heraeus, was used for the neutron scattering experiments. The Raman 
scattering measurements were performed using a standard Raman laser system. Depolarized 
spectra in 90° scattering geometry were collected in the -300-=-5000 cm -1 frequency range. 
Such a wide range was necessary to take properly into account the shape and the tempera- 
ture dependence of the luminescence background which is weak at room temperature, but 
becomes important for a quantitative determination of the coupling function at low fre- 



quency (< 20 cm -1 ) and low temperature (< 70 K). In our case the maximum luminescence 
contribution to the intensity in the problematic range < 7 cm -1 was estimated to be ~ 7% at 
12 K and w 1% at 50 K. The details of the data analysis and of the luminescence subtraction 
will be described in a more extended paper |[24|| . The neutron scattering data were taken at 
two different spectrometers (the time of flight spectrometer IN6 of the ILL at Grenoble, and 
a triple axis spectrometer at the Oak Ridge High Flux Reactor). The data underwent the 
usual corrections (subtraction of the empty container, normalization of different detectors 
to a vanadium measurement), and the density of vibrational states was derived following 
the procedure described in detail in Ref. |[25|| . Standard microcanonical molecular dynamics 
(MD) simulations of vitreous silica were performed on systems of N = 1536 and N = 5148 
atoms (box lengths L ^2.9 and 4.3 nm respectively) interacting via the two- and three-body 



potential proposed by Vashista et al. [26]. The long range interaction was treated by the ta- 



pered reaction field method [^7j] and the equations of motion were integrated by the leap-frog 



|28| alghoritm with a time step At = 0.5 fs (more details on the MD can be found in 



We used the normal mode analysis in the harmonic approximation and MD to derive the 
dynamical quantities of interest, i.e. the density of vibrational states g(oS) and the Raman 
spectrum Ir{uj). The computation of the Raman scattering needs a model for the effec- 
tive atomic polarizability and its dependence on the relative atomic position. In this study 
we associated a point polarizability to each atom (ao/cesi — 4) and used the electric dipole 
propagator (T^(r)) to describe the effective atomic polarizability, IL = Y,j oiiOCjT^ {fi — fj). 
Even though this polarization model may be a crude approximation it gives a reasonable 
agreement with the experimental data. Further details of the calculation will be given in 
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In Fig. 1 we report the Raman and neutron reduced spectra J(R t N)(w,T): 

Jw)(w,T) -,[nM+ir () 

that, according to Eq. (1), correspond to g{u)C {u) / u 2 . The data reported in Fig. 1 show 
that the frequency dependence of the spectra obtained by the two techniques is similar at 



all the investigated temperatures. Indeed, at room temperature, the maximum of the boson 
peak (uj bp ) is centered at oj bp ~ 45 cm -1 and 35 cm -1 in the RS and INS spectra respectively 
and the small blue shift of the RS maximum is due to the coupling function. Moreover, the 
temperature dependence of uj bp was recently investigated in connection with the results of 
inelastic x-ray scattering study |I2"| |. By increasing T, the BP frequencies, as seen by RS 



and INS, increase at the same rate, and at the same rate as the collective excitations at 
Q = 1.6 nm _1 observed by inelastic X-ray scattering. This temperature dependence was 
temptatively ascribed to a next-nearest neighbour fourth-order anharmonicity [|12[ . 

The ratios of Raman to neutron data collected at different temperatures in the range 
50-j-llOO K are reported in Fig. 2 where it is evident that the data at all temperatures are 
practically coincident for to >k, 30 cm -1 , while the presence of non negligible QES produces 
a marked T dependence in the low frequency region. The QES has been from time to 
time assigned to relaxation process, anharmonicity or multi-phonons scattering; its origin 
is beyond the purpose of the present work. As a consequence of QES, it is only in the 
T — > limit that the Ir/In ratio can be identified, apart from a proportionality factor, 
with the coupling function C(u). However, by reducing the temperature, also in this low 
frequency region the data appear to pile up on the continuation of the straight line that fits 
the frequency behaviour in the 30-80 cm -1 range. The inset in Fig. 1 shows that above 
« 20 cm -1 the 50 and 12 K spectra coincide, indicating that QES is negligible at these 
frequencies. 

In Fig. 3, we compare the MD simulated coupling function with the experimental one at 
the lowest temperature. The MD data have the advantage of being free from QES contribu- 
tions, but below 30 cm -1 , owing to the finite size of the simulation box, the statistics is too 
poor to give reliable results ||30|| . In the reported range the agreement between simulated and 
experimental results is very satisfactory. As mentioned, the experimental data at T = 51 
K were fitted by a straight line in the 30-80 cm -1 range (full line in Fig. 3). We note that 
our coupling function does not follow a simple power-law, and seems to extrapolate to a 
non-vanishing value at uj — > as emphasized by the dashed lines in Fig. 3 that delimit the 



±2cr confidence band of the fit. 

The non vanishing value of C(u — > 0) gives information on the characteristics of the 
vibrational modes in topologically disordered systems. It is well know that in crystalline 
materials the eigenvectors of the modes are plane waves, i. e. ej(p) oc exp(ik p ■ r}) where 
j labels the atoms and p the different modes which have well defined wavenumbers k p . In 
the acoustic region the coupling function behaves as C(u) ~ uj 2 [Q. This uo 2 law can be 
understood by considering that the scattering intensity is proportional to the square of the 
relative displacements of atoms, the main contribution coming from nearest neighbours; for 
the p-th mode this relative displacement is proportional to |e}(p) — e}_i(p)| ~ |Ve.,-(p)| oc k p , 
and hence to u> p /v. At variance with the crystalline case, in a topologically disordered system 
the eigenvectors of the modes are no longer pure plane waves. As suggested by recent MD 



calculations on a model LJ glass |pl| , |32| and on a model vitreous silica ||33|| , the eigenvector 
of the p-th mode consists of a plane wave part, with a rather well defined wavenumber k p , 
and of a further component €j, whose characteristics are not correlated with k p , and strongly 



resembling a white "noise" in the atomic displacements, e}(j>) = eyl — a p exp(ik p -fj) + a p ej. 
The plane-wave part accounts for the (rather) well defined peak observed in the dynamical 
structure factor of glasses, while the random-like part is the one mainly responsible for the 
value of C(oj). Indeed, at low frequencies, the local strain associated with the random-like 
component is much larger than that associated with the plane wave component. The same 
scenario is found in our MD simulation of V-S2O2. In Fig. 4 we report as an example the 
z component of the eigenvector of a normal mode (u> p ~ 80 cm -1 ) as a function of the z 
atomic coordinate in the N = 1536 atoms system. The plane wave component of this mode 
corresponds to a longitudinal mode, and superimposed to it there is the noisy part. Similar 
results concerning the characteristics of the normal modes in Y-S1O2 were recently obtained 
by Taraskin and Elliott (Fig. 4 of ref. |33|]). Whether the weight of this noisy component 
(cr p ) remains finite even at macroscopic wavelength, or, on the contrary, becomes negligible 
at k p — > 0, is still matter of discussion. The MD, due to the finite size of the simulation box, 
cannot access the very low frequeny region, and therefore no information can be obtained on 
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the magnitude of the noisy-to-plane wave ratio at long wavelengths. Indirect information on 
this issue can be gathered from the behaviour of C(u): in fact, a finite value for C(ou — > 0), as 
obtained in this work, supports the conjecture that the noisy component does not disappear 
even in the long wavelength limit. In fact, this can be the only origin of a non vanishing 
coupling coefficient at low frequency, because as mentioned above, the contribution of the 
plane wave component to C(u), which is proportional to u 2 , vanishes in the ui — > limit. 

Such conjecture is also supported by recent Brillouin light scattering studies of the tails 
of the inelatic peaks in the dynamic structure factors S(Q,u) of glasses. In Brillouin light 
scattering (Q « 0.02 — 0.04 nm" 1 ) p4| , |35| the S(Q, to) shows a tu-independent tail on the low 
frequency side of the inelastic peaks, whose value is Q- independent. This implies that, no 
matter how low the Q value is, the eigenvectors of all the modes at frequency smaller than 
the peak frequency, Q p = vQ, contain in their space Fourier transform such a Q component. 
In other words, the eigenvector of a mode at frequency Q, contains in its space Fourier 
transform all the Q components larger than Q p = Q/v. This is verified down to the lowest 
experimentally investigated frequency (uo ~ 0.001 cm -1 J34J) corresponding to a wavelength 
of ~ 10 fim. 

This result seems to contradict a well established, and textbook- reported, idea: "At 
wavelengths A, much higher than the interatomic distance a and the local order correlation 
length £, the sound wave feels the disordered system as homogeneous and all happens as 
in a continuous elastic medium" . The present results indicate that this picture needs to be 
revised at a microscopic level. Even when A/£ is large, the modes are not pure plane waves: 
the local disorder still implies non-regular local atomic displacements. 



We thank S. Elliott for a preprint of Ref. [33] and for useful discussions. 
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CAPTIONS 



1 - (a): Reduced Raman spectra, I (u) / {u)[n(u) + 1]}, in the - 150 cm— 1 range at different temper- 
atures (from top to bottom 45, 323, 423, 523, 873, and 1073 K). (b): Density of states divided by 
J 1 as obtained by inelastic neutron scattering in the frequency range as in a) at different tempera- 
tures (from top to bottom 51, 318, 523, 873, and 1104 K). The dashed line is the calculated Debye 
contribution as obtained using the sound velocities determined by Brillouin scattering at 50 K. (c): 
Reduced Raman spectra (in arbitrary units) at 12 and 50 K 

2 - C{uS) obtained by the Raman to Neutron scattering ratio at some significant temperatures (from 
top to bottom 51 (o), 224 (•), 271 (x), 318 (o), 423 (A), 523 (*), 873 (V), and 1104 (+) K). The 
dashed line is the linear fit to the 51 K data in the 30-80 cm -1 range. The inset shows the low 
frequency part of C{uo). 

3 - Comparison between C(w) obtained by simulation (full circles) and experiments at T=51 K (open 
circles). The full line is the linear fit to the experimental data in the 30-80 cm -1 range, and the 
dashed lines delimit the ±2(7 confidence band. 

4 - The z-componcnts of the eigenvector of the mode at w = 80 cm -1 , obtained by MD of S1O2, are 
reported versus the z-component of the atomic coordinates. The continuous line is the sine best fit. 
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